DCR Pseudo-section Simulation



In [1]:

    
from SimPEG import *
from SimPEG.Examples import DC_PseudoSection_Simulation

from ipywidgets import interactive, FloatText, FloatSlider, ToggleButtons #interactive plots!

%matplotlib inline



In [2]:

    
fig, ax = DC_PseudoSection_Simulation.run()









    



Transmitter 8 of 9 -> Time:0.747934818268 sec Transmitter 8 of 9
Forward completed






    



/Users/sgkang/anaconda/lib/python2.7/site-packages/matplotlib/collections.py:590: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  if self._edgecolors == str('face'):

Now that we have a tool to forward model data in 3D space, we can experiment with different survey configurations.

We here give the two options: pole-dipole (pdp) or dipole-dipole (dpdp).

In both cases we need to specify three important parameter.

a: Transmitter and receivers seperation (m)

b: Dipole seperation (m)

n: Number of receiver dipoles along line

We can also specify values of conductivity for the background (sig0) and the two sphere anomalies (sig1 and sig2).



In [3]:

    
Simul_Fct = lambda a,b,n, sig0, sig1, sig2, stype: DC_PseudoSection_Simulation.run(param = np.r_[a,b,n], sig=np.r_[sig0,sig1,sig2], stype = stype)

interactive(Simul_Fct, a=FloatText(min=10.,max=40.,step=5.,value=30.),
                        b=FloatText(min=10.,max=40.,step=5.,value=30.),
                        n=FloatText(min=1,max=30,step=5,value=10.),
                        sig0=FloatText(min=1e-4,max=1e+4,step=1e+1,value=1e-2),
                        sig1=FloatText(min=1e-4,max=1e+4,step=1e+1,value=1e-1),
                        sig2=FloatText(min=1e-4,max=1e+4,step=1e+1,value=1e-3),
                        stype=ToggleButtons(options=['pdp','dpdp'],value='dpdp'))









    



Transmitter 8 of 9 -> Time:0.968999862671 sec Transmitter 8 of 9
Forward completed






    












    





(<matplotlib.figure.Figure at 0x38c8400>,
 <matplotlib.axes._subplots.AxesSubplot at 0x187e9828>)



In [ ]:



In [ ]: